FERMI NATIONAL ACCELERATOR LABORATORY 



FERMILAB-TM-2403-E 
TEVEWWG/top 2008/01 
CDF Note 9225 
D0 Note 5626 

6th March 2008 



Combination of CDF and D0 Results 
on the Mass of the Top Quark 



The Tevatron Electroweak Working Groups 
for the CDF and D0 Collaborations 



Abstract 

We summarize the top-quark mass measurements from the CDF and D0 experiments at 
Fermilab. We combine published Run-I (1992-1996) measurements with the most recent 
preliminary Run-II (2001-present) measurements using up to 2.1 fb" 1 of data. Taking 
correlated uncertainties properly into account the resulting preliminary world average 
mass of the top quark is M t = 172.6 ± 0.8(stat) ± l.l(syst) GeV/c 2 , assuming Gaussian 
systematic uncertainties. Adding in quadrature yields a total uncertainty of 1.4 GeV/c 2 , 
corresponding to a relative precision of 0.8% on the top-quark mass. 



1 The Tevatron Electroweak Working Group can be contacted at tev-ewwg@fnal.gov. 
More information can be found at [EttpT77tevewwg.fnal.gov 
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1 Introduction 



The experiments CDF and D0, taking data at the Tevatron proton-antiproton collider lo- 
cated at the Fermi National Accelerator Laboratory, have made several direct experimental 
measurements of the top-quark pole mass, M t . The pioneering measurements were based on 
about 100 pb" 1 of Run-I (1992-1996) data [Il[2l[3l[llEl[6l[Tl[8l[9l UHl [HI Ll2l and include 
results from the it — *> qq'bqq'b (all-j), the it —> Ivqq'bb (1+j), and the it — > l+vbt'Vb (di- 
1) decay channels^]. The Run-II measurements summarized here are the most recent results 
in the 1+j, di-1, and all-j channels using 1.9 — 2.1 ffcT 1 of data and improved analysis tech- 
niques [ia[T2[l5l[T6l[I71[l8l[T9l[20l[2T]. 

This note reports the world average top-quark mass obtained by combining five published 
Run-I measurements [21 El El LTl UQl EH] with one published Run-II CDF result [13], three 
preliminary Run-II CDF results [131 [T5l [T6] and three preliminary Run-II D0 results [TTJ [181 
|2~T| . The combination takes into account the statistical and systematic uncertainties and their 
correlations using the method of references [221 [23] and supersedes previous combinations [2H 
[251 [261 [271 [28] . 

The input measurements and error categories used in the combination are detailed in Sec- 
tion [2] and [31 respectively. The correlations used in the combination are discussed in Section H] 
and the resulting world average top-quark mass is given in Section [5j A summary and outlook 
are presented in Section [61 



2 Input Measurements 

For this combination twelve measurements of M t are used, five published Run-I results, and 
one published plus six preliminary Run-II results, all reported in Table [H In general, the Run-I 
measurements all have relatively large statistical uncertainties and their systematic uncertainty 
is dominated by the total jet energy scale (JES) uncertainty. In Run-II both CDF and D0 take 
advantage of the larger it samples available and employ new analysis techniques to reduce both 
these uncertainties. In particular the JES is constrained using an in-situ calibration based on 
the invariant mass of W — > qq' decays in the 1+j and all-j channels. The Run-II D0 analysis in 
the 1+j channel constrains the response of light-quark jets using the in-situ W — ► qq' decays. 
Residual JES uncertainties associated with 7]— and pr-dependencies as well as uncertainties 
specific to the response of 6-jets are treated separately. Similarly, the Run-II CDF analysis 
in the 1+j and all-j channels also constrain the JES using the in-situ W — > qq' decays. Small 
residual JES uncertainties arising from rj— and pr-dependencies and the modeling of 6-jets are 

2 Here £ = e or /i. Decay channels with explicit tau lepton identification are presently under study and are 
not yet used for measurements of the top-quark mass. 
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Run-I published 


Run-II preliminary 






CDF 




D0 






CDF 






D0 






all-j 


1+j 


di-1 


1+j 


di-1 


1+j 


di-1 


all-j 


bey 


1+j /a 


1+j/b 


di-1 


Result 


186.0 


176.1 


167.4 


180.1 168.4 


172.7 


171.2 


177.0 


180.7 


170.5 


173.0 


173.7 


iJES 


0.0 


0.0 


0.0 


0.0 


0.0 


1.3 


0.0 


1.8 


0.0 


0.0 


0.0 


0.0 


aJES 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.7 


0.8 


1.9 


bJES 


0.6 


0.6 


0.8 


0.7 


0.7 


0.4 


0.1 


0.1 


0.0 


0.2 


0.1 


0.9 


cJES 


3.0 


2.7 


2.6 


2.0 


2.0 


0.5 


1.7 


0.6 


0.0 


0.0 


0.0 


2.1 


dJES 


0.3 


0.7 


0.6 


0.0 


0.0 


0.1 


0.1 


0.1 


0.0 


1.7 


1.4 


0.9 


rJES 


4.0 


3.4 


2.7 


2.5 


1.1 


0.2 


1.8 


0.5 


0.3 


0.0 


0.0 


0.0 


Signal 


1.8 


2.6 


2.8 


1.1 


1.8 


0.6 


0.7 


0.6 


1.4 


1.0 


0.5 


0.8 


BG 


1.7 


1.3 


0.3 


1.0 


1.1 


0.6 


0.4 


1.0 


7.2 


0.5 


0.4 


0.6 


Fit 


0.6 


0.0 


0.7 


0.6 


1.1 


0.2 


0.6 


0.6 


4.2 


0.1 


0.2 


0.9 


MC 


0.8 


0.1 


0.6 


0.0 


0.0 


0.4 


0.7 


0.3 


0.7 


0.0 


0.0 


0.2 


UN/MI 


0.0 


0.0 


0.0 


1.3 


1.3 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


Syst. 


5.7 


5.3 


4.9 


3.9 


3.6 


1.7 


2.8 


2.4 


8.5 


2.2 


1.7 


3.4 


Stat. 


10.0 


5.1 


10.3 


3.6 


12.3 


1.2 


2.7 


3.3 


14.5 


1.9 


1.3 


5.4 


Total 


11.5 


7.3 


11.4 


5.3 


12.8 


2.1 


3.9 


4.1 


16.8 


2.9 


2.2 


6.4 



Table 1: Summary of the measurements used to determine the world average M t . All numbers 
are in GeV/c 2 . The error categories and their correlations are described in the text. The 
total systematic uncertainty and the total uncertainty are obtained by adding the relevant 
contributions in quadrature. 

included in separate error categories. The Run-II CDF di-1 measurement uses a JES determined 
from external calibration samples. Some parts of the associated uncertainty are correlated with 
the Run-I JES uncertainty as noted below. 

In previous combinations the Run-II CDF 1+j analysis used the JES determined from the 
external calibration as an additional Gaussian constraint. This required us to treat that mea- 
surement as two separate inputs in the combination in order to accurately account for all the 
JES correlations. This Gaussian constraint is not used in the present analysis as it does not 
significantly improve the sensitivity. Thus we can treat this measurement as a single input in 
the same manner as all the other measurements. 

As discussed in the previous combination, a new analysis technique from CDF is included 
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(lxy). This measurement uses the mean decay- length from B-tagged jets to determine the 
top-quark mass. While the statistical sensitivity is not nearly as good as the more traditional 
methods, this technique has the advantage that since it uses only tracking information, it is 
almost entirely independent of JES uncertainties. As the statitistics of this sample continue to 
grow, this method could offer a nice cross-check of the top-quark mass that's largely indepen- 
dent of the dominant JES systematic uncertainty which plagues the other measurements. The 
statistical correlation between this measurement and the Run-II CDF 1+j measurement is deter- 
mined using Monte Carlo signal-plus-background psuedo-experiments which correctly account 
for the sample overlap and is found to be consistent with zero (to within < 1%) independent 
of the assumed top-quark mass. 

The two D0 Run-II lepton+jets results [T7J [TS] are derived from Run-IIa and Run-IIb 
datasets, respectively, and are labelled as such. The D0 Run-II dilepton result [2T] is itself a 
combination of two results using different techniques but the same di-lepton data set (191 [20] . 

Table [1] also lists the uncertainties of the results, sub-divided into the categories described 
in the next Section. The correlations between the inputs are described in Section HJ 



3 Error Categories 

We employ the same error categories as used for the previous world average [28]. They include 
a detailed breakdown of the various sources of uncertainty and aim to lump together sources of 
systematic uncertainty that share the same or similar origin. For example, the "Signal" category 
discussed below includes the uncertainties from ISR, FSR, and PDF - all of which affect the 
modeling of the it signal. Additional categories are included in order to accommodate specific 
types of correlations. For example, the jet energy scale (JES) uncertainty is sub-divided into 
several components in order to more accurately accommodate our best estimate of the relevant 
correlations. Each error category is discussed below. 

Statistical: The statistical uncertainty associated with the M t determination. 

iJES: That part of the JES uncertainty which originates from in-situ calibration procedures 
and is uncorrelated among the measurements. In the combination reported here it cor- 
responds to the statistical uncertainty associated with the JES determination using the 
W — > qq' invariant mass in the CDF Run-II 1+j and all-h measurements. Residual JES 
uncertainties, which arise from effects not considered in the in-situ calibration, are in- 
cluded in other categories. 

a JES: That part of the JES uncertainty which originates from differences in detector e/h 
response between fe-jets and light- quark jets. It is specific to the D0 Run-II measurements 
and is taken to be uncorrelated with the D0 Run-I and CDF measurements. 
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bJES: That part of the JES uncertainty which originates from uncertainties specific to the 
modeling of 6-jets and which is correlated across all measurements. For both CDF and 
D0 this includes uncertainties arising from variations in the semi-leptonic branching 
fraction, 6-fragmentation modeling, and differences in the color flow between 6-jets and 
light-quark jets. These were determined from Run-II studies but back-propagated to the 
Run-I measurements, whose rJES uncertainties (see below) were then corrected in order 
to keep the total JES uncertainty constant. 

cJES: That part of the JES uncertainty which originates from modeling uncertainties corre- 
lated across all measurements. Specifically it includes the modeling uncertainties associ- 
ated with light-quark fragmentation and out-of-cone corrections. 

dJES: That part of the JES uncertainty which originates from limitations in the calibration 
data samples used and which is correlated between measurements within the same data- 
taking period, such as Run I, Run Ha or Run lib, but not between experiments. For CDF 
this corresponds to uncertainties associated with the r/-dependent JES corrections which 
are estimated using di-jet data events. For D0 Run-II this corresponds to uncertainties 
associated with the light-quark response as determined using the W —>■ qq' invariant mass 
in the 1+j channel and propagated to the di-1 channel. The residual ^-dependent and 
Pr-dependent uncertainties for the D0 Run-II measurements are also included here since 
they are constrained using Run-II 7+jet data samples. 

r JES: The remaining part of the JES uncertainty which is correlated between all measurements 
of the same experiment independent of data-taking period, but is uncorrelated between 
experiments. This is dominated by uncertainties in the calorimeter response to light- 
quark jets. For CDF this also includes small uncertainties associated with the multiple 
interaction and underlying event corrections. 

Signal: The systematic uncertainty arising from uncertainties in the modeling of the it signal 
which is correlated across all measurements. This includes uncertainties from variations 
in the ISR, FSR, and PDF descriptions used to generate the it Monte Carlo samples 
that calibrate each method. It also includes small uncertainties associated with biases 
associated with the identification of 6-jets. 

Background: The systematic uncertainty arising from uncertainties in modeling the domi- 
nant background sources and correlated across all measurements in the same channel. 
These include uncertainties on the background composition and shape. In particular 
uncertainties associated with the modeling of the QCD multi-jet background (all-j and 
1+j), uncertainties associated with the modeling of the Drell-Yan background (di-1), and 
uncertainties associated with variations of the fragmentation scale used to model W+jets 
background (all channels) are included. 

Fit: The systematic uncertainty arising from any source specific to a particular fit method, 
including the finite Monte Carlo statistics available to calibrate each method. 
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Monte Carlo: The systematic uncertainty associated with variations of the physics model 
used to calibrate the fit methods and correlated across all measurements. For CDF it 
includes variations observed when substituting PYTHIA [29J [301 EI] (Run I and Run II) 
or ISAJET [32J (Run I) for HERWIG [331 EI] when modeling the it signal. Similar 
variations are included for the D0 Run-I measurements. The D0 Run-II measurements 
use ALPGEN [35] to model the it signal and the variations considered are included in 
the Signal category above. 

UN/MI: This is specific to D0 and includes the uncertainty arising from uranium noise in the 
D0 calorimeter and from the multiple interaction corrections to the JES. For D0 Run-I 
these uncertainties were sizable, while for Run-II owing to the shorter integration time 
and in-situ JES determination, these uncertainties are negligible. 



These categories represent the current preliminary understanding of the various sources of 
uncertainty and their correlations. We expect these to evolve as we continue to probe each 
method's sensitivity to the various systematic sources with ever improving precision. Varia- 
tions in the assignment of uncertainties to the error categories, in the back-propagation of the 
bJES uncertainties to Run-I measurements, in the approximations made to symmetrize the 
uncertainties used in the combination, and in the assumed magnitude of the correlations all 
negligibly effect (<g 0.1GeV/c 2 ) the combined M t and total uncertainty. 



4 Correlations 



The following correlations are used when making the combination: 



• The uncertainties in the Statistical, Fit, and iJES categories are taken to be uncorrelated 
among the measurements. 

• The uncertainties in the aJES and dJES categories are taken to be 100% correlated among 
all Run-I and all Run-II measurements on the same experiment, but uncorrelated between 
Run I and Run II and uncorrelated between the experiments. 

• The uncertainties in the rJES and UN/MI categories are taken to be 100% correlated 
among all measurements on the same experiment. 

• The uncertainties in the Background category are taken to be 100% correlated among all 
measurements in the same channel. 

• The uncertainties in the bJES, cJES, Signal, and Generator categories are taken to be 
100% correlated among all measurements. 
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Run-I published 


Run-II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j lxy 


D0 

1+j/a 1+j/b di-1 


CDF-I 1+j 
CDF-I di-1 
CDF-I aU-j 


1.00 

0.29 1.00 

0.32 0.19 1.00 








D0-I 1+j 
D0-I di-1 


0.26 0.15 0.14 
0.11 0.08 0.07 


1.00 

0.16 1.00 






CDF-II 1+j 
CDF-II di-1 
CDF-II aU-j 
CDF-II lxy 


0.30 0.17 0.16 
0.45 0.27 0.33 
0.17 0.11 0.15 
0.11 0.03 0.02 


0.22 0.09 
0.21 0.11 
0.09 0.05 
0.10 0.01 


1.00 

0.24 1.00 

0.11 0.17 1.00 

0.16 0.03 0.02 1.00 




D0-II 1+j/a 
D0-II 1+j/b 
D0-II di-1 


0.16 0.09 0.06 
0.11 0.06 0.04 
0.18 0.12 0.11 


0.11 0.05 
0.09 0.03 
0.17 0.08 


0.16 0.07 0.06 0.10 
0.12 0.04 0.04 0.09 
0.14 0.19 0.07 0.01 


1.00 

0.20 1.00 

0.21 0.14 1.00 



Table 2: The resulting matrix of total correlation coefficients used to determined the world 
average top quark mass. 

Using the inputs from Table [1] and the correlations specified here, the resulting matrix of total 
correlation co-efficients is given in Table [2j 

The measurements are combined using a program implementing a numerical x 2 minimiza- 
tion as well as the analytic BLUE method [22l|23]. The two methods used are mathematically 
equivalent, and are also equivalent to the method used in an older combination (36], and give 
identical results for the combination. In addition, the BLUE method yields the decomposition 
of the error on the average in terms of the error categories specified for the input measure- 
ments [23] • 



5 Results 



The combined value for the top-quark mass is: 

M t = 172.6 ± 1.4 GeV/c 2 , (1) 

with a x 2 °f 6.9 for 11 degrees of freedom, which corresponds to a probability of 81% indicating 
good agreement among all the input measurements. The total uncertainty can be sub-divided 
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Run-I published 


Run-II preliminary 


CDF 
1+j di-1 all-j 


D0 
1+j di-1 


CDF 

1+j di-1 all-j lxy 


D0 

1+j /a 1+j/b di-1 


Pull 

Weight [%] 


+0.5 -0.5 +1.2 
-4.2 -0.7 -0.6 


+1.5 -0.3 
+1.8 +0.2 


+0.1 -0.4 +1.2 +0.5 
+35.8 +9.7 +8.8 -0.7 


-0.8 +0.2 +0.2 
+ 15.2 +35.2 -0.6 



Table 3: The pull and weight for each of the inputs used to determine the world average mass 
of the top quark. See Reference [22] for a discussion of negative weights. 



into the contributions from the various error categories as: Statistical (±0.8), total JES (±0.9), 
Signal (±0.5), Background (±0.4), Fit (±0.1), Monte Carlo (±0.2), and UN/MI (±0.02), for a 
total Systematic (±1.1), where all numbers are in units of GeV/c 2 . The pull and weight for each 
of the inputs are listed in Table [31 The input measurements and the resulting world average 
mass of the top quark are summarized in Figure [D 



The weights of many of the Run-I measurements are negative. In general, this situation 
can occur if the correlation between two measurements is larger than the ratio of their total 
uncertainties. This is indeed the case here. In these instances the less precise measurement 
will usually acquire a negative weight. While a weight of zero means that a particular input 
is effectively ignored in the combination, a negative weight means that it affects the resulting 
central value and helps reduce the total uncertainty. See reference [22] for further discussion of 
negative weights. 

Although the x 2 from the combination of all measurements indicates that there is good 
agreement among them, and no input has an anomalously large pull, it is still interesting to 
also fit for the top-quark mass in the all-j, 1+j, and di-1 channels separately. We use the same 
methodology, inputs, error categories, and correlations as described above, but fit for the three 
physical observables, M t all ~ J , Ml + \ and M t dl_1 . The results of this combination are shown in 
Table [Hand have x 2 °f 4.2 for 9 degrees of freedom, which corresponds to a probability of 90%. 
These results differ from a naive combination, where only the measurements in a given channel 
contribute to the M t determination in that channel, since the combination here fully accounts 
for all correlations, including those which cross-correlate the different channels. Using the 
results of Tabled] we calculate the chi-squared consistency between any two channels, including 
all correlations, as x 2 (dil — Ij) = 0.8, x 2 i}j ~ odlj) — 1-5, and x 2 ( a ^j ~ dil) = 2.7. These 
correspond to chi-squared probabilities of 39%, 23%, and 10%, respectively, and indicate that 
the determinations of M t from the three channels are consistent with one another. 
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Parameter 


Value (GeV/c 2 ) 


Correlations 


M all-j 

M t 1+j 
M? 1 - 1 


177.3 ±3.9 

172.4 ±1.5 
169.8 ±3.1 


1.00 

0.12 1.00 

0.18 0.26 1.00 



Table 4: Summary of the combination of the nine measurements by CDF and D0 in terms of 
three physical quantities, the mass of the top quark in the all-jets, lepton+jets, and di-lepton 
channel. 

6 Summary 

A preliminary combination of measurements of the mass of the top quark from the Tevatron 
experiments CDF and D0 is presented. The combination includes five published Run-I mea- 
surements and one published plus six preliminary Run-II measurements. Taking into account 
the statistical and systematic uncertainties and their correlations, the preliminary world-average 
result is: M t = 172.6 ± 1.4 GeV/c 2 , where the total uncertainty is obtained assuming Gaussian 
systematic uncertainties and adding them plus the statistical uncertainty in quadrature. While 
the central value is somewhat higher than our 2007 average, the averages are compatible as 
appreciably more luminosity and refined analysis techniques are now used. 

The mass of the top quark is now known with a relative precision of 0.8%, limited by 
the systematic uncertainties, which are dominated by the jet energy scale uncertainty. This 
systematic is expected to improve as larger data sets are collected since new analysis techniques 
constrain the jet energy scale using in-situ W — >• qq' decays. It can be reasonably expected 
that with the full Run-II data set the top-quark mass will be known to much better than 
0.8%. To reach this level of precision further work is required to determine more accurately 
the various correlations present, and to understand more precisely the 6-jet modeling, Signal, 
and Background uncertainties which may limit the sensitivity at larger data sets. Limitations 
of the Monte Carlo generators used to calibrate each fit method may also become important as 
the precision reaches the ~ 1 GeV/c 2 level and will warrant further study in the near future. 
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Mass of the Top Quark (*Preliminary) 



Measurement 


OUr-l 


ai-i 


r\rx I 
\JVj-\ 


Hi 1 

ai-i 


nnr ii 


*41 1* 

ai-i 


r-\rx 1 1 

D0-II 


^i: i* 

di-l 


CDF-I 


l+j 


D0-I 


l+j 


CDF-II 


l+j* 


D0-II 


l+j/a* 


D0-II 


l+j/b* 


CDF-I 


all-j 


CDF-II 


all-j* 


CDF-II 


Ixy 



Tevatron Run- 




M top [GeV/cl 
167.4 ± 11.4 
1 68.4 ± 12.8 
171.2 ± 3.9 
173.7± 6.4 
176.1 ± 7.3 
180.1 ± 5.3 
172.4± 2.1 
170.5± 2.9 
173.0± 2.2 

- 186.0 ±11.5 
177.0 ± 4.1 

- 180.7 ±16.8 



% I dot = 6.9 / 1 1 
172.6± 1.4 



150 170 

M top [GeV/c 2 ] 



190 



March 2008 



Figure 1: A summary of the input measurements and resulting world average mass of the top 
quark. 
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